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ABSTRACT 


An  equation  is  derived  for  the  growth  in  time  of  a  barium 
cloud  in  a  strongly  conducting  background.  The  expression  is 
correct  to  first  order  in  an  expansion  in  powers  of  the  ratio  of 
the  cross -cloud  conductivity  to  the  background  conductivity. 
Numerical  solutions  of  this  equation  indicate  significant  dif¬ 
ferences  from  the  zero  order  result.  The  forward  moving  part 
of  the  cloud  thins  and  indents  while  the  backside  steepens. 

The  list  of  weakly  ionized  instabilities  is  extended  to 
include  the  true  gravitational  instability.  This  is  the  closest 
contender  to  the  E  X  B  but  still  appears  to  have  a  threshold 
higher  by  about  a  factor  of  ten. 
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I .  INTRODUCTION 

1  2 

In  the  previous  semi-annual  technical  reports  *  we  have 
discussed  the  various  instabilities  which  a  low'/3  ,  weakly- 
ionized  plasma  is  subject  to  and  demonstrated  that  the  threshold 
for  onset  of  the  E  X  B  instability  is  lower  than  for  all  others. 
In  order  to  further  identify  the  instability,  it  is  necessary  to 
develop  a  realistic  model  of  the  equilibrium  state  of  the  cloud 
just  before  onset  of  striation.  A  model  sufficient  for  this 
purpose  was  developed  and  reported  in  I,  Linearization  of  the 
equations  of  motion  about  this  equilibrium  state  then  provides  a 
convenient  set  of  equations  for  numerical  calculation.  Some  of 
the  results  obtained  were  reported  in  II.  Further  results  are 
noted  in  Section  IV  of  this  report.  We  are  now  designing  an 
improved  numerical  scheme  which  will  allow  considerable  increase 
in  grid  size  and  should  advance  our  knowledge  of  the  critical 
eigenvalues  and  eigenstates. 

In  the  meantime  emphasis  has  shitted,  for  the  moment,  to 
extending  the  model  of  the  equilibrium  state.  The  previous  model 
assumed  a  background  ionosphere  which  was  infinitely  conducting. 
This  meant  that  there  was  no  dielectric  reduction  of  the  ambient 
electric  field  in  the  interior  of  the  cloud  due  to  charges 
collecting  on  the  surface.  As  a  result  the  cloud  moved  smoothly 
(gaussian  density  distribution)  and  with  a  velocity  equal  to  that 
of  the  ions  alone  (modified  only  by  the  internal  ambipolar 
electric  field) .  We  have  now  extended  this  model  to  include  the 
first  order  correction  in  an  expansion  in  powers  of  the  small 
ratio  of  the  cross-cloud  conductivity  to  the  ionosphere’s 
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conductivity.  This  extension  allows  dielectric  screening  of 
the  external  field  to  occur.  Time -dependent  numerical  solutions 
of  these  modified  equilibrium  equations  give  evidence  of  flat¬ 
tening  of  the  rear  of  the  cloud  and  indentation  and  thinning  in 
the  direction  of  drift. 

The  extended  equilibrium  model  is  discussed  in  Section  II 
below.  Numerical  results  are  given  in  Section  III.  A  brief 
discussion  of  the  instability  eigenvalue  calculations  is  in 
Section  IV  as  well  as  extension  of  previous  instability  esti¬ 
mates  to  include  the  "true"  gravitational  instability. 
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II.  EQUILIBRIUM  STATE  IN  A  HIGHLY-CONDUCTING  BACKGROUND 

In  the  rest  frame  of  the  neutrals,  the  equation  of 
motion  of  the  barium  plasma  takes  the  form 

c'W/  -  Di  Vi-W-  Dh  V7W  ±  ui  51-  CWE^ 

'SIC  y 

±  Jdn  7,  C M  El.')  +-  b  •  y  w  X  E  =  O 

(1) 

where  the  symbols  all  have  the  same  meaning  as  in  I.  if 

one  now  subtracts  the  electron  form  of  Equation  (1)  from  the  ion 

form,  we  obtain 

- ( ol- or ) v7  w  4  co,7-D,n  1-  (^r4^)&-(tv§) 

+  (JU»  +jtU)  V„  (WE,,)  —  ( (strT/xI-  (A'tVjJil)  b- 

(2) 

In  the  equilibrium  state  we  assume  a  smooth  profile  with 
characteristic  length jv  and  X  in  the  radial  and  parallel  direction 

of  the  order  of  ^  Di  j  ol\  and  with  Ex  =.  —  $  /ft  t  Ec^ 

H\i  ~  —  4^/jL  where  is  the  ambient  electric  field 

in  the  neutral  rest  frame.  We  also  assume  Jl  (j>  ^  fcT  where 
T  denotes  the  common  temperature.  Dividing  all  terms  by  Du  /$/" 
and  using  the  Einstein  relations  and  the  usual  in¬ 
equalities  (5UT\Z  (Sir)*  9  JXw  ^  /An  »  Di  ^  Oi. 

it  is  easy  to  see  that  all  terms  in  Equation  (2)  except  for  the 
second  and  fourth  are  of  order  Du/Du  I  •  We 

assume  C  Eo  5-  lyC1'  0 *  order  ten  or  lsss  and  t^us  not 
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large  enough  to  change  the  ordering  due  to  VA,;  which  is 
approximately  (1/1300).  We  also  assume/?  and  X  of  similar  order 
and  {/7t)+c£  order  unity  or  larger. 

To  lowest  order  in  an  expansion  in  powers  of  /&\[ 
our  equation  becomes: 

Vii  [  A  V  +  ~  0 

where 

^  =  Ai //%  ~  KT~/(L  (3) 

Integration  yields 

A  V„N  .  (4) 

Assuming  symmetry  of  the  equilibrium  about  the  s?  - o  plane,  we 
have  ^  ~  6  •  Integrating  once  more  yields 

A  a/  ~  f  +  3  (%A  ~  °  (S) 

Unlike  the  method  used  in  I ,  we  shall  confine  our  expansion  in 
K  /  3,7  to  the  lowest  order,  so  that  Equation  (5)  is 

our  basic  equation  for  determining  E.  However,  again  unlike  I, 
we  will  not  assume  an  infinitely  conducting  background  but  will 
use  the  following  improved  method  to  determine  £  • 

By  our  general  assumption  of  quasineutrality,  we  must 

have 

£■•  £  =  °  (6) 

at  every  point  in  space,  where  d  denotes  the  electrical  current. 
If  we  now  integrate  Equation  (6)  along  a  magnetic  field  line 
passing  through  the  cloud,  from  a  point  infinitely  far  from  the 
cloud  on  one  side  to  a  similar  point  on  the  other  side,  and  use 
the  vanishing  of  the  parallel  flux  at  infinity,  one  obtains 
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(7) 


f  A-jTj.  <3  *  ~  f  £.•£ 

<r  tcoj)  '-'to/vo 

Here  we  have  separated  the  surviving  terms  into  a  portion  inte¬ 
grated  over  the  cloud  and  that  integrated  over  the  ionosphere. 
By  Equation  (S) ,  we  see  that  the  perpendicular  electric  field 
is  independent  of  j?  outside  the  cloud  (Note  that  ^  H  is  con¬ 
stant  outside  the  cloud).  As  a  result,  the  r.h.s.  of  Equation 
(7)  may  be  written  as 


-L 


aP  Zj_  ■  £4.  * 


9- 


It  *46 


or 


-  2T. 


B 


P,  •  £ 


where  is  the  Pedersen  conductivity  in  the  ionosphere  and 
is  its  height  integral. 

We  thus  obtain  the  equation 


The  integral  on  the  r.h.s.  is  proportional  to  0^  ,  the  charac¬ 
teristic  cross  field  conductivity  of  the  cloud,  and  thus  the 
right  hand  side  is  proportional  to  / Zq  where  <tis  a 

characteristic  width  of  the  cloud.  We  assume  that  the  cross 
field  conductivity  in  region  B  is  very  large  compared  to  the 
cloud  conductivity,  hence  this  parameter  is  small.  If  one  now 
expands  E j_  in  powers  of  \  jone  immediately  obtains 

the  zero-order  solution  _  f'O 

Zl  ■  §J-  =  0  ■ 

Since  the  external  field  must  be  equal  to  the  ambient  field 

(o)  __  p- 

far  from  the  cloud,  this  yields  E.j_  ~  Ez  c  . 

We  now  use  this  result  to  calculate  £  in  the  cloud, 
and  then  N  in  the  cloud.  By  Equation  (5) 
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-  A  Mk  N  f-  j{&x ) 


E  -  -  A  Vti 
/ V 


t- 


E 


(9) 


where  we  assume  the  first  term  vanishes  outside  the  cloud.  Now 
substitute  Equation  (9)  in  the  fourth  and  sixth  terms  of  Equa¬ 
tion  (1)  and  then  eliminate  the  fifth  term  by  the  appropriate 
combination  of  the  ion  and  electron  forms  of  Equation  (1) .  The 
result  is: 

-  cD„  KM  -  *Dt_  &V  /  Cl- «0  PV 


3  t 


<-  c-§-  (k,  V)  M  - 

-O 


O 


(10) 


This  equation  is  identical  to  Equation  (12)  of  I  if  one  trans¬ 
forms  to  the  frame  moving  with  velocity  ccc  -  All 

symbols  have  the  same  meaning  as  in  this  equation.  Thus  the 
solution  for  N  follows  immediately  from  Equation  (13)  of  X  and  is 

v  EK P  jLL  ? 

I  ^IT'  *  "  J  .  (11) 

assuming  an  impulse  source  of  strength  ^  at  the  origin  at 
The  fact  that  the  resulting  density  distribution  is 
identical  with  that  in  I  is  not  surprising  since  to  lowest  order 
in  A  »  the  background  ionosphere  appears  to  be  infinitely  con¬ 
ducting  (relative  to  the  cross  cloud  conductivity).  It  is  also 
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no  surprise  that  the  cloud  is  cylindrically  symmetric  since 
there  can  be  no  dielectric  distortion  in  this  same  approxima¬ 
tion.  (The  only  charge  that  can  build  up  on  the  outer  surface 
of  the  cloud  is  one  that  is  uniform  over  the  entire  surface.) 
This  is  the  nature  of  the  ambipolar  field  specified  by  the 
first  term  in  Equation  (9).  Note  that  J^is  a  constant  along 
the  field  lines,  as  assumed. 

Let  us  now  proceed  to  first  order  in  A  .  Equation  (8) 


becomes 


The  zero-order  current  is  obtained  at  once  by  using  the  usual 
diffusion  expression  together  with  Eqs.  (9)  and  (11).  The 
result  takes  a  particularly  simple  form  if  one  makes  the  for¬ 
mal  transformation 

A?  =  +  &~-f 

^vv  ***** 

t'  -  t  . 


)  Bat  <JC,K  ,4  t 


(12) 


In  terms  of  these  coordinates  the  expression  above  takes  the 
form 


where  we  have  omitted  the  primes,  set  6^  -  -  p^^and  assumed 
&  .  The  boundary  condition  is  that  vanish  far 

from  the  cloud. 
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Using  the  Green's  function  method,  the  solution  is: 

</>  W,  O  ■  ztF  ^  ^  O 

where  /t*  =  and  ^  denotes  the  r.h.s.  of  Equation  (13). 

Using  the  well-known  expansion  (see  Morse  and  Feshbach,  Equa¬ 
tion  (10.1.18)  ^  a vf 

s€y,  £  =  Ms  ^  "  Z.  ™  ^r~)  [M  ( t-?')] 

where  is  the  greater  of  A±  and  Jl'x  and  is  the  lesser,  it  is 
easy  to  carry  out  the  angular  integral.  The  result  is: 


r  - 


-  * 

rr  ZB 


n- 


t- 


*KP(j£f)  *'u'- 

One  can  now  calculate  most  of  the  remaining  integrals  with  the 


result 


roO 

r 


f 


j 

A 


-f 


B 

c 


•AT  [  &rj+ 
r  £*PC  v3£t) J 


A. 
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The  resulting  electric  field  follows  by  differentiation.  We 
obtain 


(14) 
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where  C-n  is  a  unit  vector  in  the  radial  direction.  The  obvious 
extension  of  Equation  (9)  is 


£ 


0) 


(15) 


If  we  also  transform  Equation  (1)  to  the  coordinates  expressed 
by  Equation  (12) ,  substitute  Equation  (15)  in  the  fourth  and 
sixth  terms  and  eliminate  the  fifth  term  between  the  electron 
and  ion  forms  of  this  equation,  the  result  is 


3AL  -  J>„  1 7„\y  -<Pj_  ^4-/9  Vl-  tNE'n) 

Sir  & 

■h  £■  O-r)  k  ■  VA/y  E,n=  O 

(16) 

(l) 

where  £  is  given  by  Equation  (14) .  This  is  the  equation  of 
motion  of  the  ion  density  as  seen  from  a  frame  moving  relative 
to  the  neutrals  with  the  transformation  given  by  Equation  (12), 
i.e...a  frame  moving  with  the  velocity  of  the  cloud  in  an  in¬ 
finitely  conducting  background. 
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III.  NUMERICAL  SOLUTION  FOR  TIME  DEVELOPMENT  OF  THE 
EQUILIBRIUM  STATE 

Equation  (16)  is  reduced  to  two-dimensional  form  by 
averaging  over  the  z -direction.  Thus  Equation  (16)  reduces  to 

-  JTj.  Vx^  +  p,  ( N  |l" ) 

[i7) 

where  [Q-  .  A  computer  solution  is  then  obtained 

by  assuming  an  initial  solution 

....  -**/«.’"• 
n(o)=  jk_  e. 

Tra1- 

measuring  all  lengths  in  terms  of  a  and  time  in  units  of  Cl 
(The  proper  value  of  jzCl'  is  then  given  by  Equation  (14)  with  % 
replaced  by  Cu/i •)  Tw0  basic  dimensionless  parameters 

also  appear  in  the  dimensionless  form  of  Equation  (17).  These 
are: 

^  ~ 

If  a1 

which  is  the  ratio  of  the  cross-cloud  conductivity  to  that  of  the 
background  and 

Y-  CL  Eo<l 

fctTr+T-) 

which  is  the  normalized  external  electric  field  strength.  The 
solution  is  obtained  by  a  difference  scheme  using  the  A.D.I. 
method.  Further  details  will  be  given  in  a  separate  report. 

Typical  results  are  shown  in  Figures  1-5.  Here  one 

chose 

CL  -  H-.  S’  Km 

"Dj.  ~  Ckm) 
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and  the  length  scale  is  then  shown  in  units  of  kilometers. 

The  magnetic  field  is  out  of  the  paper  and  JE0  is  in  the  posi¬ 
tive  y-  direction.  Note  the  interesting  distortion  which  occurs 
in  time.  The  rear  of  the  cloud  tends  to  flatten  and  steepen, 
while  the  front  side  indents  and  thins.  All  these  results  have 
a  simple  interpretation  in  terms  of  the  polarization  produced 
in  the  cloud  by  Pedersen  currents  (as  well  as  Kail  currents  at 
low  (fl^) . 
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IV. 


OTHER  WORK  IN  PROGRESS 


The  numerical  scheme  for  evaluating  the  stability 
boundaries  and  for  determining  the  shape  of  the  eigenfunctions 
is  being  completely  revised.  The  previous  method,  described  in 
II,  was  pushed  to  an  11x11  grid  size  with  no  significant  change 
in  the  general  form  of  the  results.  Now  we  plan  to  use  a  much 
faster  method  which  yields  only  the  most  rapidly  growing  eigen¬ 
value,  instead  of  all  eigenvalues.  Advice  and  assistance  is 
being  obtained  from  the  Los  Alamos  Scientific  Laboratory. 

The  catalog  of  low  frequency  instabilities  in  weakly 
ionized  plasma  was  examined  in  I  to  determine  the  most  likely 
candidate  for  the  origin  of  barium  cloud  striation.  The  E  X  B 
was  shown  to  have  a  thr  shold  about  a  factor  of  100  lower  than 
that  for  the  ion  sound  instability  (either  E  X  B  driven  or  driven 
by  diamagnetic  streaming)  or  the  drift  dissipative  instability. 

The  gravitational  instability  due  to  magnetic  field  curvature  is 

negligible  compared  to  the  above;  however  the  ’’true”  gravitational 

* 

instability  is  a  serious  contender.  We  have  investigated  the 
competition  between  this  mode  and  the  E  X  J3  by  solving  the 
stability  problems  in  similar  idealized  slab  geometries  (Kt=0  ). 
Roughly  speaking,  the  E  X  B  goes  unstable  when 

Eo  >  _Du_ 

while  the  gravitational  is  unstable  for 


/mt  ^  ^ 


>UTi 


*  We  are  indebted  to  Dr.  Francis  Perkins  for  this  observation. 
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Here  L  is  roughly  o*  '].*  order  of  the  characteristic  length  cf 
the  density  gradient  (although  it  tends  to  be  somewhat  smaller 
for  the  gravitational  case).  Comparing  the  two  results  above, 
we  see  that  the  E  X  B  has  the  lower  threshold  if 

U.  =  ce.  > H  +  (•5ir)t] 

B  Sl+ 

Since  JV+  =  24  sec  ,  'j.  «  9.8m/sec  and  Uo  is  of  the  order  of 
20-100  m/sec  in  most  releases,  we  see  that  the  E  X  B  dominates  by 
about  a  factor  of  10  if  (jtT')+=,2.  as  observations  indicate. 
Clearly  this  instability  is  the  closest  contender.  It  must  be 
remembered,  however,  that  there  have  been  releases  in  which  the 
field  lines  were  nearly  vertical  and  which  still  demonstrated 
rapid  formation  of  striations. 
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